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ABSTRACT 

We reexamine likelihood analyzes of the Local Group (LG) acceleration, paying particular 
attention to nonlinear effects. Under the approximation that the joint distribution of the LG 
acceleration and velocity is Gaussian, two quantities describing nonlinear effects enter these 
analyzes. The first one is the coherence function, i.e. the cross-correlation coefficient of the 
Fourier modes of gravity and velocity fields. The second one is the ratio of velocity power 
spectrum to gravity power spectrum. To date, in all analyzes of the LG acceleration the second 
quantity was not accounted for. Extending our previous work, we study both the coherence 
function and the ratio of the power spectra. With the aid of numerical simulations we obtain 
expressions for the two as functions of wavevector and o-g ■ Adopting WMAP's best determina- 
tion of a^, we estimate the most likely value of the parameter (3 and its errors. As the observed 
values of the LG velocity and gravity, we adopt respectively a CMB-based estimate of the LG 
velocity, and Schmoldt et al.'s (1999) estimate of the LG acceleration from the PSCz catalog. 
We obtain [3 = 0.661^;^^; thus our errorbars are significantly smaller than those of Schmoldt 
et al. This is not surprising, because the coherence function they used greatly overestimates 
actual decoherence between nonlinear gravity and velocity. 

Key words: methods: numerical, methods: analytical, cosmology: theory, dark matter, large- 
scale structure of Universe 



1 INTRODUCTION 

Analyzes of large-scale structure of the Universe provide esti- 
mates of cosmological parameters that are complementary to those 
from the cosmic microwave background (CMB) measurements. 
In particular, comparing the large-scale distribution of galaxies 
to their peculiar velocities enables one to constrain the quantity 
P = /b. Here, f2,„ is the cosmological matter density param- 
eter and b is the linear bias of galaxies that are used to trace the 
underlying mass distribution. This is so because the peculiar veloc- 
ity field, V, is induced gravitationally and therefore is tightly cou- 
pled to the matter distribution. In the linear regime, this relationship 
takes the form 

- = n°jJ^Sir)^, (1) 

where 5 denotes the mass density contrast and distances have been 
expressed in km ■ s~^. Under the assumption of linear bias S = 
b~^Sg, where Sg denotes the density contrast of galaxies, and the 
amplitude of peculiar velocities depends linearly on f3. 

These comparisons are done by extracting the density field 
from full-sky redshift surveys (such as the PSCz; Saunders et 
al. 2000), and comparing it to the observed velocity field from pe- 
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culiar velocity surveys. The methods for doing this fall into two 
broad categories. One can use equation Q to calculate the pre- 
dicted velocity field from a redshift survey, and compare the result 
with the measured peculiar velocity field; this is referred to as a 
velocity-velocity comparison. Alternatively, one can use the dif- 
ferential form of this equation, and calculate the divergence of the 
observed velocity field to compare directly with the density field 
from a redshift survey; this is called a density-density comparison. 

Peculiar velocities of galaxies and groups of galaxies are gen- 
erally determined by measuring their distances independently of 
redshifts. However, the motion of the Local Group (LG) of galax- 
ies can be deduced in another way, namely from the observed 
dipole anisotropy of the CMB temperature. This dipole reflects, 
via the Doppler shift, the motion of the Earth with respect to the 
CMB rest frame. The components of this motion of local, non- 
cosmological origin (the Earth motion around the Sun, the Sun 
motion in the Milky Way and the motion of the Milky Way in the 
LG) are known and can be subtracted (e.g., Courteau & van den 
Bergh 1999). When transformed to the barycenter of the LG, the 
motion is towards {l,b) — (276° ± 3° , 30° ±2°), and of amplitude 
■"LG ~ 627 ± 22 km • as inferred from the 4-year COBE data 
(Lineweaver et al. 1996). It can be compared to that predicted from 
a redshift survey and historically this was the first velocity-velocity 
comparison (with angular data only, e.g. Meiksin & Davis 1986, 
Yahil, Walker & Rowan-Robinson 1986; with redshift information, 
Strauss & Davis 1988, Lynden-Bell, Lahav & Burstein 1989). This 
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velocity estimate is much more accurate tlian estimates of peculiar 
velocities based on redshift-independent distances. Therefore, the 
analysis of the LG motion remains an interesting alternative to cur- 
rent velocity-velocity comparisons, performed simultaneously for 
many galaxies, but with less accurately measured velocities. 
Let us define the scaled gravity, g, by the equation: 



g 



d r r 



(2) 



where again distances have been expressed in km • s^^. The scaled 
gravity is proportional to the gravitational acceleration, and can be 
measured from a redshift survey. Equation Q yields 



(3) 



so a velocity-velocity comparison is in fact a velocity-gravity one. 
Hereafter, we will refer to 'scaled gravity' simply by 'gravity'. It is 
also often called 'the clustering dipole'. 

A number of effects (nonlinear effects, observational windows 
through which the velocity and gravity of the LG are observed, 
shot noise) spoil the linear relationship 0. Therefore, to estimate 
P one cannot simply equate the two dipoles, but a more sophisti- 
cated approach is needed. Commonly adopted is a maximum like- 
lihood approach, which provides a way to account for these effects 
and to compute errors of estimated cosmological parameters. Here 
we reexamine a likelihood analysis of Schmoldt et al. (1999; S99), 
paying particular attention to proper modelling of nonlinear effects 
(NE), that affect the estimated values of parameters and their errors. 

In a previous work (Chodorowski & Ciecieljg 2002; hereafter 
C02) we concentrated on one quantity describing NE in the LG 
velocity-gravity comparison, the coherence function (CF). It is the 
cross-correlation coefficient of the Fourier modes of the gravity and 
velocity fields. We showed that the form of the function adopted by 
S99 drastically overestimates actual decoherence between nonlin- 
ear gravity and velocity. This implies that the true random error of 
P is significantly smaller. In the present study we give a complete 
description of NE present in the analysis. In particular, we show 
that not only the CF is relevant, but also the ratio of the power 
spectrum of velocity to the power spectrum of gravity. With aid of 
numerical simulations we model the two quantities as functions of 
the wavevector and of cosmological parameters. We then combine 
these results with observational estimates of vlg and gLG, and 
obtain the 'best' value of (3 and its errors. The paper is organized 
as follows. In Section|2|we outline an analytical description of the 
likelihood for cosmological parameters, based on the measured val- 
ues of the LG velocity and acceleration. In Section |3| we describe 
our numerical simulations, and we model numerically the coher- 
ence function and the ratio of the power spectra. An estimation of 
the most likely value of /3 and its errors is presented in Section|4] 
Summary and conclusions are in Section|5| 



2 ANALYTICAL DESCRIPTION OF THE LIKELIHOOD 

Let /(g, v) denote the joint distribution function for the LG grav- 
ity and velocity. It is commonly approximated by a multivariate 
Gaussian (Strauss et al. 1992, hereafter S92; S99). This approxima- 
tion has support from numerical simulations (Kofman et al. 1994, 
Ciecieljg et al. 2003), where the measured nongaussianity of g and 
V is small. This is rather natural to expect since, e.g., gravity is an 
integral of density over effectively a large volume, so the central 
limit theorem can at least partly be applicable (but see Catelan & 
Moscardini 1994). 



In a Bayesian approach, one ascribes a priori equal probabili- 
ties to values of unknown parameters, which allows one to express 
their likelihood function via /: 



/I(param.) = f(v,g | param. 



(4) 



As the parameters to be estimated we adopt /3 and b. Using statisti- 
cal isotropy of g and v, the distribution function can be cast to the 
form (Juszkiewicz et al. 1990, Lahav, Kaiser & Hoffman 1990): 



/(g,v) = 



(1 



-3/2 



■ exp 



2rfixy 



2(1 



(5) 



where Ug and ctv are the r.m.s. values of a single Cartesian com- 
ponent of gravity and velocity, respectively. From isotropy, Ug = 
(g • g) /3, and crj = {v • v)/3, where {■) denote the ensemble aver- 
aging. Next, (x, y) = (g/o-g, v/ctv), and /i — cosi/; with?/) being 
the misalignment angle between g and v. Finally, r is the cross- 
correlation coefficient of g,n with Vm, where g,n (vm ) denotes an 
arbitrary Cartesian component of g (v). From isotropy. 



Also from isotropy. 



(6) 



(7) 



where denotes the Kronecker delta. In other words, there are 
no cross-correlations between different spatial components. 

We have to take into account that the LG gravity is measured 
through the window Wg : 



If V is irrotational like g, then similarly 

4-K 



(8) 



(9) 



where 9 is the (minus) velocity divergence, 8 = —V ■ v, and Wv is 
the velocity window. Due to Kelvin's circulation theorem, the cos- 
mic velocity field is vorticity-free as long as there is no shell cross- 
ing. N-body simulations (Bertshinger & Dekel 1989, Mancinelli 
et al. 1994, Pichon & Bernardeau 1999) show that the vorticity of 
velocity is small in comparison to its divergence even in the fully 
nonlinear regime. 

The best current estimate of the LG gravity is inferred from 
the PSCz catalog of galaxies (Saunders et al. 2000). S99 fol- 
low S92 and measure the LG gravity through the standard IRAS 
window. 



(r/rsf 



W^^i 1, 
0, 



r < rs , 
Vs <r < R 

max ; 

R max 



(10) 



The window is characterized by a small-scale smoothing and a 
sharp large-scale cutoff. Following S99, we adopt the values rs = 
500 km ■ s~^ and -Rmax = 15, 000 km ■ s^^, appropriate for the 
PSCz catalog. The window function relevant to the LG velocity is 







1, Otherwise , 



(11) 



which has a small-scale cutoff, rmin ~ 100 km ■ s~^, to reflect 
the finite size of the LG (S92). Using numerical simulations, we 
have checked that the velocity field remains approximately Gaus- 
sian even for such a small smoothing scale. 

In Fourier space, relations Js} and j9} read: 



3 



(12) 



(13) 



where the subscript k denotes the Fourier transform. The quantity 
W is not a Fourier transform of W, but is related to it in the fol- 
lowing way (S92): 



W{k) = kl W(r)ji{kr)dr. 



(14) 



Here and below ji represents the spherical Bessel function of order 
I. In particular, 

W,{k) = -jo(fcfl^,,), (15) 

krs 

and 

VK,(fc) = jo(fcr^in). (16) 
From equations <12t and <13t we have 
1 1"^ ^ 

and 

{v-v)=2^y Wl{k)Pe{k)dk . (18) 

Here, P(fc) and P0{k) are respectively the power spectrum of the 
density and the power spectrum of the velocity divergence. Thus, 



(v-v) 



where 



7^(fc) 



27r2 



W^{k)TZ{k)P{k)dk, 



P.{k) _ Pe(k) 



(19) 



(20) 



Pg(fc) P(k) 

and Pv and Pg are the power spectra respectively of velocity and 
gravity. Furthermore, 



v) 



1 

2^ 



M7g(fc)W^v(fc)C(fc)Pg/^(A:)P'/^(fc)dfc, (21) 



where C{k) is the so-called coherence function^ (S92), or the cor- 
relation coefficient of the Fourier components of the gravity and 
velocity fields: 

(gk • vi^) _ {5^61) 



^^^^ ~ (|gkP)V2(|vk|2)l/2 (j5^|2)l/2^|^^|2)l/2 

Hence, we obtain 

/o°° Ws{k)w^{k)c{k)n^''^{k)P{k)dk 



1/2 



/o°° Wiik)P{k)dk Wd{k)TZ{k)P{k)dk 



1/2 



(22) 



.(23) 



Equations <17> . <19t and J23> specify all parameters (the vari- 
ances and the correlation coefficient) that determine distribution j5} 
in the absence of observational errors. The deviation of the correla- 
tion coefficient from unity is then due to different windows, through 
which the gravity and the velocity of the LG are measured, and due 
to nonlinear effects. The latter are described by two functions: the 

^ S92 call it the decoherence function. We prefer the name 'coherence', 
because higher values of the function imply higher, not lower, correlation 
between gravity and velocity. 



coherence function, and the ratio of the power spectra. In the linear 
regime, C(fc) = 1 and TZ = . This yields 



Wg(fc)Wv(fc)P(fc)dfc 



1/2 



J^Wi{k)P{k)dk W^{k)P{k)dk 



1/2 



(24) 



so for linear fields the correlation coefficient is determined solely 
by the windows, as expected. 



3 MODELLING NONLINEAR EFFECTS 
3.1 Numerical simulations 

We follow the evolution of the dark matter distribution using the 
pressureless hydrodynamic code CPPA (Cosmological Pressure- 
less Parabolic Advection, see Kudlicki, Plewa & Rozyczka 1996, 
Kudlicki et al. 2000 for details). It employs an Eulerian scheme 
with third order accuracy in space and second order in time, which 
assures low numerical diffusion and an accurate treatment of high 
density contrasts. Standard applications of hydrodynamic codes in- 
volve a coUisional fluid; however, we implemented a simple flux in- 
terchange procedure to mimic collisionless fluid behaviour. Thanks 
to this approach we avoid a few problems of N-body codes. The 
main advantage of a hydrodynamic code over an N-body code 
is accurate treatment of low density regions. Moreover it directly 
produces a volume-weighted velocity field. This is important be- 
cause in the definition of the CF, equation j22> . the velocity field 
is volume-weighted, not mass-weighted. Furthermore, the Fast 
Fourier Transform can be directly applied to the data on an uni- 
form grid. 

The computational domain forms a (400 /!.~^Mpc)'^ cube 
with 256'^ grid cells and periodic boundary conditions. This setup 
allows us to cover a broad range of wavenumbers, k £ [0.016, 2.] 
/i/Mpc, which is important for accurate calculation of the integrals 
in equation <23> . The initial distribution of the density fluctuations 
is Gaussian and their power spectrum is given by a CDM formula 
(as in Eq. 7 of Efstathiou, Bond & White 1992) with the shape pa- 
rameter r — 0.19, as inferred from the IRAS PSCz survey (Suther- 
land et al. 1999) and in agreement with recent determinations from 
WMAP (Spergel et al. 2003). The initial value of the square root of 
the variance of the density field in spheres of radius 8 Mpc, 
crg(ti), is 0.019. We studied a range of outputs from the simu- 
lations, corresponding to different values of as- Specifically, full 
output data were dumped in constant intervals of the scale factor, 
Aa = 0.025. 

All our simulations assume Sim ~ 0.3, f^A ~ 0. In the 
mildly non-linear regime, the quantity v = n~'' ''v is insensitive to 
the cosmological density parameter and cosmological constant, as 
demonstrated both analytically (Bouchet et al. 1995, Nusser & Col- 
berg 1998; see also Appendix B3 of Scoccimarro et al. 1998) and 
by means of N-body simulations (Bernardeau et al. 1999). There- 
fore, our results should be valid for any cosmology. Specifically, 
we have 



7^ = Qli^TZ 



and 

C = C, 
where 



TZik) 



Pv(fc) 



(25) 



(26) 



(27) 
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and 



Cik) 



|gkl2)l/2(|vkl2)l/2 



(28) 



In a previous paper (C02), we tested numerically the dependence of 
the coherence function on Q,m and found it to be extremely weak. 

In C02, we also investigated numerical effects. In short, we 
found that the effects of resolution affect numerical determination 
of the CF at scales smaller than 4 grid cells, ie. twice the Nyquist 
wavelength. Therefore, all the results we present here are for < 1 
/i/Mpc. For these wavenumbers, both simulated C and TZ practi- 
cally do not depend on resolution. Also, here we use a larger box 
size than in C02 in order to model longer modes, while the short 
wavelength limit (k < 1 /i/Mpc) still allows us to accurately cal- 
culate the integrals in Equation <23> . as will be shown in subsec- 
tionlsH 



3.2 Coherence function 

We studied the CF in C02. We found there that the characteris- 
tic decoherence scale is an order of magnitude smaller than previ- 
ously used (S92). The weak point of the formula fitted in C02 was 
a poor description of the CF for long-wavelength modes. These 
modes are however important when calculating integrals in equa- 
tion <23> . (From eQs. ll2ltT3l it follows that the gravity and velocity 
fields are more sensitive to long-wavelength modes than the den- 
sity field.) Therefore, here we use another fitting function, which is 
more accurate for low values of k: 



C{k) = [l + (fflofc - aafc^-^ + a^k^ f 



(29) 



Parameters a; were obtained for 35 different values of erg in the 
range [0.1,1], and we found the following, power-law, scaling rela- 
tions: 



ao 
tti 
a2 



4.908 al'"^° 

o ceo 0-734 

2.663 (Jg 
5.889 (7g 



(30) 



The fit was calculated for k G [0, 1] /i/Mpc, with the imposed 
constraint C(fc = 0) = 1. Unlike the previous fit, the new one 
results in the value of the correlation coefficient of gravity and ve- 
locity which agrees with that measured directly in our simulations. 
In FigureQwe show, for various values of erg, the CF from simula- 
tions, its fit <29t , and results of perturbative calculations, described 
in C02. We see that the perturbative approximation breaks down for 
erg > 0.5. 



3.3 Ratio of tlie power spectra 

The ratio of the velocity to the gravity power spectra, TZ, is related 
to its scaled counterpart, TZ, by equation <25> . The quantity TZ, de- 
fined in equation <27> . practically does not depend on the back- 
ground cosmological model. It departs from unity in the nonlinear 
regime because the velocity grows slower than it would be expected 
from the linear approximation. 

We have found that TZ obtained from simulation can be fitted 
with the following formula: 



7^(fc) = [1 + (7.071/0)-*]-°, 
with 

a = -0.06574 + 0.29195crg 



for 0.3 < ag < 1 



(31) 



(32) 



• •»imjjUj^;_j^|^ 1 |..., 


-J—J — i — ^ — — 1 '■■ 1 

o 


1 

=0.2 






=0.51 






=0.87^> 


1,1.1,1 


1,1.1,1 


1 , 



0.3 0.4 0.5 0.6 
k [h/Mpc] 



Figure 1. The coherence function for various values of erg . Dots show the 
results of simulations, solid lines - our fit <29l . and dashed lines - the results 
of perturbative calculations. Perturbative approximation breaks down for 
o-g > 0.5. 
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Figure 2. The ratio of the velocity power spectrum to the gravity power 
spectrum as a function of the wavevector for three different values of erg . 
Dotted lines show the results of our simulations and solid lines - the fit I31i . 



We stress that the above relation between a and ag is valid for 
erg G [0.3, 1], which is still a sufficiently wide range of values. 
Figure|2|shows the ratio of the power spectra from simulation, and 
our fit. 



3.4 Convergence of r 

We calculate the correlation coefficient, r, by inserting fits <29> and 
<31> into formula <23> . However, the integrals in this formula ex- 
tend over the whole fc-space, while the fits have been obtained for a 
limited range of wavenumbers between 0.016 and 1 /i/Mpc. There- 
fore, the question has to be answered whether this extrapolation is 
justified. 

Figure|3| shows integrands of the integrals in the formula <23> 
for (Tg = 0.84. It is evident that contributions from wavenumbers 
greater than unity are negligible. This is so because observational 
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Figure 3. Integrands of the integrals appearing in equation 1231 as func- 
tions of the wavevector for erg = 0.84. Since the fc-axis is logarithmic, 
the ordinate is multiplied by and extra power of fc, so equal areas under a 
function correspond to equal contributions to an integral. Solid line shows 
the product WgW^CTZ^^^Pk, dotted - W^UPk, and dashed - W^Pk. 
Units of the ordinate axis are arbitrary. 



cosmological component, (7g^, and errors, e . Since gravity here is 
inferred from a galaxian, rather than mass, density field, we have 
(jg J, = where Sg is the variance of a single component of the 
true (i.e., mass-induced) gravity. From equation <17> . 



1 



Wg{k)P{k)dk. 



(33) 



The gravity errors are twofold: due to finite sampling of the 
galaxy density field, and due to the reconstruction of the galaxy 
density field in real space. Therefore, = ((jIjv, + a'^^c) /3, where 
cr|p^ and (T^cc are respectively the shot noise (or, sampling) variance 
and the reconstruction variance. (Both a"^^ and a'^^^ are full, i.e., 
3D, variances.) Estimated by S99 using mock catalogs, the cumu- 
lative shot noise at 150 h^^ Mpc amounts to (jsn = 160 km ■ s^^. 
An average reconstruction error in the differential contribution to 
the cumulative gravity, produced by a shell of matter IQ h^^ Mpc 
wide, is 15 km • s~^. Since up to 150 Mpc there are 15 such 
shells, we have (Jroc = \/T5 ■ 15 km • 58 km • s~^. To sum 

up. 



2 2 

,22 I f^SN + ""rcc 
Sg -t- g 



where 



o"sN ~ 160 km ■ s 



and 



= 58 km • s 



(34) 



(35) 



Errors in the measured velocity of the LG are negligible com- 
pared to those in the gravity. Equations <19t and <25> yield 



— ^^n 



(36) 



windows of the LG gravity and velocity filter out smaller scales.'^ 
This is not quite the case for wavenumbers smaller than 0.016 
/i/Mpc, but these scales are well within the linear regime, for which 
the limiting values of the CF and the ratio of the power spectra are 
known to converge to unity. 



4 PARAMETER ESTIMATION 

We are now applying our formalism to the PSCz survey. As the 
value of (78 we adopt its WMAP's estimate, erg — 0.84 (±0.04; 
Spergel et al. 2003). This specifies the coherence function and the 
ratio of the power spectrum of velocity to the power spectrum of 
gravity. Also, this provides a normalization for the power spectrum 
of density. 



4.1 Parameter dependence of the model 

As stated before, the likelihood of specific values of /3 and b is de- 
termined by the distribution js}- In this distribution, the observables 
are g and v, or g, v, and the misalignment angle, tp. Following 
S99, we adopt for them the following values: g — 933 km ■ s^^ 
(from the distribution of the PSCz galaxies up to 150 Mpc), 
V = 627 km • s^^ (inferred from the 4-year COBE data by 
Lineweaver et al. 1996), and tl) = 15°. 

The theoretical quantities are cTg, CTv, and r. The variance of a 
single spatial component of measured gravity, (jg, is a sum of the 



^ Strictly speaking, the velocity window passes contributions from 
wavenumbers up to about 2 /i/Mpc, but the ratio of the power spectra 
damps them additionally for fc > 1 /i/Mpc. 



where 



1 

6^ 



Wi(k)1Z(k)P{k)dk. 



(37) 



Finally, errors in the estimate of the LG gravity do not affect the 
cross-correlation between the LG gravity and velocity, but increase 
the gravity variance. This has the effect of lowering the value of the 
cross-correlation coefficient. Specifically, from equations <25> - J26> 
and <34t we have 



2 2 

362 s| 



-1/2 



where 



j^W^{k)W^{k)C{k)'R}/'^{k)P{k)dk 



1/2 r 



/o°° wi{k)P{k)dk wd(k)n{k)P{k)dk 



1/2 



(38) 



.(39) 



Thus, the likelihood depends explicitly on the parameters b and 
ilm. However, for reasons that will become evident later on, as 
the parameters to be estimated we choose b and /3 = Q^lf /b. Since 



— (5b, both CTg, CTv and r depend on b. On the other hand. 



only (Jv depends on (5. This makes the derivation of the most likely 
value of 13, given b, simple. 



4.2 /3 for known b 

Using equations ^5), <36> . and the equality 



■ /3b, the logarith- 



mic likelihood for (3 takes on the form: 

ln£(/3) = -31n27r - 31n[ag6sv(l - r^)^/^] - 31n/3 



1 f g^ Ir^gv 



2(1 -r2) 



(3'^b'^si ag(3bs^ 



(40) 
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Figure 4. Likelihood contours for tfie pai'ameters /3 and b, con'esponding 
to the confidence levels of 68, 90 and 95%. The maximum of the likelihood 
function is denoted by a dot. 



To find its maximum, we calculate its partial derivative with respect 
to (3, dlnC/dp, and equate it to zero. This yields the following 
equation: 



3(1 



= 0. 



(41) 



The LG gravity, inferred from the PSCz survey, is tightly coupled 
to its velocity, 1 — r ^ 1 and 1 — ^ <c; 1. (Specifically, /i = 0.97 
and, for erg around 0.8, r ~ 0.93.) At first approximation we can 
therefore assume r = /i = 1, hence 



Pi — 7 . 

osv g 

Using equation <34t we obtain finally 

1/2 



/3i = 



1 + 



2 2 



(42) 



(43) 



Thus, the best estimate of j3 is not just the ratio of the LG velocity 
to its gravity: it is modified by nonlinear effects (which affect Sv 
through the function TZ), different observational windows (which 
affect differently Sg and Sv), and observational errors. 

At next approximation, in equation <41> one could approxi- 
mate by PI . However, we have considered the case of known 
b for illustrative purposes only and from now on we relax this as- 
sumption. In the next subsection we will analyze the joint likeli- 
hood for 13 and b. 

4.3 Joint likelihood for /3 and b 

Figure|4|shows isocontours of the joint likelihood for (3 and h, cor- 
responding to the confidence levels of 68, 90 and 95%. The maxi- 
mum of the likelihood is denoted by a dot. The corresponding val- 
ues of (3 and b are respectively 0.66 and 0.94. 

This figure illustrates the well-known problem of fim-bias 
degeneracy in cosmic density-velocity comparisons. A 'common 
wisdom' is that in these comparisons, only a degenerate combi- 
nation of Q,m and b, (3 = /b, can be determined. Here this 
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Figure 5. The marginal distribution for /3. The result of marginalizing over 
all possible values of b (from zero to infinity) is shown as a dashed line. The 
result of marginalizing over the values of b in the range [0.7, 1.1] is shown 
as a solid line. 



is not strictly true, since we have estimated the most likely values 
of both f3 and b, so in principle we could solve for the best value 
of Qrn alone. However, the isocontours of the likelihood are much 
more elongated along the &-axis than along the /9-axis, making the 
resulting constraints on fim much weaker than on l3. In practice, 
therefore, from our analysis we cannot say much about ilm. or bias 
separately. On the other hand, we can put fairly tight constraints on 
/3. 

To do this, one approach is not to use any external information 
on bias. In such a case, we marginalize over all possible values of 
b (from zero to infinity). The resulting distribution for (3 is shown 
in Figure|5]as a dashed line. The distribution is positively skewed, 
with lower limit on /3 being stronger than the upper one. Specifi- 
cally, we find that [3 = 0.64l°;J* (68% confidence hmits). 

A better approach is to use the fact that the square root of the 
variance of the PSCz galaxy counts at 8 ft"^ Mpc is af^'^^ ~ 0.75 
(Sutherland et al. 1999). Combined with WMAP's estimate of as, 
this yields for the bias of the PSCz galaxies the value about 0.9.^ 
Therefore, we adopt here a conservative prior for the bias, namely 
that it is constrained to lie in the range [0.7, 1.1]. These limits are 
marked in Figure |4| as dashed vertical lines. We marginalize the 
likelihood over the values of b in this range. The resulting distri- 
bution for 13 is shown in Figure|5|as a solid line. This distribution 
is also skewed and more peaked than the previous one. We obtain 
/3 = 0.66lo:o7 (68% confidence limits). 



5 SUMMARY AND CONCLUSIONS 

We have performed a likelihood analysis of the LG acceleration, 
paying particular care to nonlinear effects. We have adopted a 
widely accepted assumption that the joint distribution of the LG ac- 
celeration and velocity is Gaussian. Then, two quantities describing 
nonlinear effects are relevant. The first one is the coherence func- 
tion, or the cross-correlation coefficient of the Fourier components 
of the gravity and velocity fields. The second one is the ratio of the 
power spectrum of the velocity to the power spectrum of the grav- 
ity. Extending our previous work, we have studied both the coher- 



Our best value of the bias, 0.94, is, given the errors, surprisingly close to 
this estimate. 
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ence function and the ratio of the power spectra. Using numerical 
simulations we have performed fits to the two as functions of the 
wavevector and . Then, we have estimated the best values of the 
parameter (3 and its errors. We have obtained (3 — 0.661q qt at 
68% confidence level. 

The analysis of the LG acceleration performed by S92 and S99 
was in a sense more sophisticated than ours. Both teams analyzed a 
differential growth of the gravity dipole in subsequent shells around 
the LG. Instead, here we used just one measurement of the total (in- 
tegrated) gravity within a radius of 150 Mpc. Nevertheless, the 
errors on /3 we have obtained are significantly smaller than those of 
S92 and S99. In particular, S99 obtained /3 = 0.70t[! Jg at Icr 
confidence level. Comparing these errors to ours should be done 
with caution, because S99 considered the joint likelihood for /3 and 
the index of the power spectrum, (To obtain a constraint on 
/3, they marginalized the distribution over the values of V allowed 
by the constraints on the Hubble constant.) Still, it is striking that 
while our best value of /3 is close to theirs, our errors are signif- 
icantly smaller. The reason is our careful modelling of nonlinear 
effects. In a previous paper (C02) we showed that the coherence 
function used by S99 greatly overestimates actual decoherence be- 
tween nonlinear gravity and velocity. Tighter correlation between 
the LG gravity and velocity should result in a smaller random error 
of /3; in the present work we have shown this to be indeed the case. 

The second factor in the LG acceleration analysis, describing 
nonlinear effects, is the ratio of the power spectra of velocity to 
gravity. Unlike the coherence function, this quantity has not been 
accounted for previously. It affects the value of the correlation coef- 
ficient, hence the random error, to lesser extent than the CF. '' How- 
ever, it does affect the most likely value of /3, as can be easily no- 
ticed from illustrative equation <42t (via Sv). Neglecting different 
nonlinear growth rates of the gravity and the velocity is equivalent 
to setting TZ = 1; we have checked that then the best value of 13 is 
0.62. Therefore, a small discrepancy between our and S99's most 
likely value of (3 is not due to nonlinear effects. 

Unlike ours, the analysis of S99 included also other con- 
straints on the velocity field around the LG: its (small) shear, and 
the value of the bulk flow within 30 h^^ Mpc. This may have 
had an effect on the best value of (3. Moreover, an approach with 
multiple windows should allow to further tighten the errors. It is 
therefore interesting to repeat the analysis of S99 exactly, but with 
proper treatment of nonlinear effects. We plan to do this in the fu- 
ture. 
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The constrains on V obtained by S99 were extremely weak, so we de- 
cided to fix the value of F and instead to study the dependence on b. 
^ If we write 7?.(A;) = 1 — e(fc), from equations I38l - I39j it is straightfor- 
ward to show that r = r.^_-^ + 0(t'^). 



